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MAGNETIZATION AND SUSCEPTIBILITY OF THE HEISENBERG FERROMAGNET 

by Lawrence Flax 
Lewis Research Center 

SUMMARY 

The anisotropic Heisenberg ferromagnet developed by Flax and Raich is reevaluated 
to include an applied magnetic field. The present calculation is restricted to the case 
of an isotropic exchange interaction, and the random phase approximation is used. The 
present theory predicts a broadening of the transition region with increasing magnetic 
field. Thermodynamic quantities, such as magnetization and susceptibility, are studied 
near the Curie point. 


INTRODUCTION 

According to the Weiss theory a ferromagnetic body possesses an internal field 
which is proportional to the magnetization. As the temperature of a ferromagnet is 
increased, the magnetization decreases until, at a temperature known as the Curie tem- 
perature T , the sample becomes paramagnetic. This change from ferromagnetic to 
paramagnetic is referred to as the ferromagnetic phase transition. Experiments and 
theory have shown that electrical, mechanical, and many thermodynamic properties of 
a material are altered when the material undergoes a phase transition. 

One way of describing the ferromagnetic phase transition is to look upon it as an 
order-disorder transition in the material. The property of order in a system describes 
the correlation between the properties of two atoms, even when the atoms are widely 
separated in space. The spontaneous magnetization present in a ferromagnet is due to 
the correlation between the spins of the atoms. Above the Curie temperature the long- 
range order is gone if there is no applied field and the system is paramagnetic. 

At absolute zero the spontaneous magnetization <S Z ) is equal to its maximum value, 
while at any other temperature T with no applied field the value depends on the ratio 
T/T . (Symbols are defined in the appendix. ) 

The effect of the external magnetic field is twofold: (1) the magnetization is in- 
creased somewhat above its zero-field value, and (2) the transition region is broadened. 



Even though temperature tends to destroy this alinement, the field causes some order- 
ing to be present. Thus, instead of an abrupt disorder at the critical temperature there 
is a tail due to some ordering. 

This report investigates some thermodynamic properties of the Heisenberg ferro- 
magnet is an external magnetic field. The objective is to obtain expressions for some 
of the thermodynamic quantities that are needed to evaluate the use of ferromagnetism 
for a refrigeration system (refs. 1 and 2). The thermodynamic quantities studied here 
are the magnetization and susceptibility. The statistical method used to determine the 
thermodynamics of the system is the double-time Green's function. Generalized Wat- 
son sums arise in the theory of magnetism and other disorder phenomena. These sums 
are evaluated analytically in order to calculate the thermodynamic quantities as a func- 
tion of temperature. Previous evaluations of these sums have been based on series ex- 
pansions and are valid only over limited temperature regions - usually away from the 
critical temperature. The methods used in references 3 to 7 and in this work enable 
one to obtain analytical expressions for the thermodynamic quantities without using 
series expansions. 


HEISENBERG FERROMAGNET 

The Heisenberg model assumes that the magnetic electrons are in states localized 
about the lattice sites with an exchange interaction taking place between the electrons. 

The model does not take into account an itinerant- electron picture. The model is 
thought to be very good for insulating ferromagnets such as europium oxide and poor for 
conductors such as iron and nickel. However, it appears that this model gives better 
results for conductors at low temperatures than those calculated from a band theory 
approach as shown by Argyle, Charap, and Pugh (ref. 8). The success of the Heisen- 
berg model is somewhat surprising since it does not take into account the spreading of 
the electronic energy levels into bands. 

The Heisenberg ferromagnet with spin 1/2 was analyzed by Bogolyubov and Tyablikov 
(ref. 9) using the techniques of double-time temperature -dependent Green functions. A 
convenient review of Green functions and Tyablikov's application of them to ferromag- 
netism is given by D. N. Zubarev in reference 10 and in references cited there. 

The Hamiltonian for the Heisenberg model is 

H - -^b h o E - E j « [<f«f * \ { s t s i + s , :s i)] <« 


where 


2 



r 


s± = s* ±i sy 


( 2 ) 


/i B is the Bohr magneton, g is the Lande g factor, Hq is the applied magnetic field 
which is assumed to be in the -z direction, is the spin operator for a spin at site i, 
is the exchange interaction between spins on sites i and j, and the sum is carried 
over all sites in the crystal. The exchange interaction is assumed to be a function of the 
distance between sites only. The self-exchange terms such as and vanish. 

Since the thermodynamic quantities of interest are the magnetization and suscept- 
ibility of the system, one is interested in evaluating correlation functions of the form 
(S^nSg) and hence of the Green's function ((S*; S^)) . 

Using the Hamiltonian (1) together with the known spin commutation relations 



one can write the equation of motion (see ref. 3 or 4) in the form 


(3) 


(E - „ B gH 0 K < S m » = i < Sp « gmA * 2 £ J gf < < s+sf; S^> ) 


lists' 


'g’ ®m' ^ 


(4) 


For a system of interacting spins, the second and third terms on the right side of equa- 
tion (4) contain Green's function with three operators and are called higher order 
Green's functions. The equation of motion for one of these higher order Green's func- 
tions involves the next higher order Green's function. In this way an infinite hierarchy 
of equations is generated. Thus, the exact treatment of the equations of motion involves 
the solution of an infinite set of coupled equations for an infinite number of Green's func- 
tions. To obtain a closed solution for a Green's function, this hierachy of equations 
must be truncated at some point. The procedure usually adopted is to make a so-called 
decoupling approximation in which the higher order Green's functions on the right side of 
equation (4) are expressed in terms of lower order Green's functions. By doing this, 
an explicit approximate solution of equation (4) for ( ( S^; S^) ) can be found. 
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RANDOM PHASE APPROXIMATION 


The essential approximation of the random phase approximation (RPA) method of 
Bogolyubov and Tyablikov (ref. 9), Tahir-Kheli and ter Haar (ref. 11), Flax and Raich 
(refs. 3 to 7), and those cited in references 3 to 7 consists of ignoring the correlations 
between (S z ) and the other operators and and replacing the operator by 
the statistical average ( S z ) , that is, 

<< S g S f’ S m>>-^ < sZ >« s J‘; S m>) (5) 

We may note that ( S z ) = (S z ) is independent of the site g on account of the translational 
invariance of the lattice. This method of truncating the higher order Green functions is 
certainly an ad hoc one since S - (S ) is not necessarily small. At low temperatures, 
one may justify the use of this decoupling procedure since the spins are nearly com- 
pletely alined. In the nearly alined state, the deviation of the spin (z component) at any 
one site from the average value would be small, and thus for a first approximation one 
might neglect - (S z ). However, the low -temperature expansion of the resulting ex- 

& 3 

pression for the magnetization contains a T term in conflict with the careful spin wave 
analysis of Dyson. 

At high temperatures the spins are more disordered, and the RPA approximation 
leads to further discrepancies: 

(1) The magnetization obeys a square root law in temperature in the vicinity of the 
Curie point as opposed to the cube root law behavior suggested by experiment. 

(2) The specific heat is zero above T c because of the abrupt disappearance of long- 
range order to contrast to experimental observation in which a nonzero value persists 
because of local or short-range order. 

However, the random phase approximation shares these defects with all theories 
based on Landau's treatment of second-order phase transitions. 

Upon applying the RPA to the equation of motion (4), one obtains 

[E - ^ B gH 0 - 2<S z )J(0)] Agm =i <^)5 gm - 2<£?> 2 J gf A fm (6) 

7/ p 1 


where for brevity 


gm 


= «S+;S m » 


(7) 
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Translational invariance allows the use of a Fourier transformation to the reciprocal 
lattice; that is, 


A m = ~ 2 exp[ik • (g - m)] A fc 
gm N k K 


( 8 ) 


and 


J(k) = J2 J gf exp [ik • (f - g)] (9) 

where N is the total number of sites in the lattice and the sums are restricted to the 
first Brillouin zone. It was shown in references 3 to 5 and 7 that the Green's function 
for this model can be written in the form 


A - < sZ ) 
k jr[E-E(k)] 


where the elementary excitation energy E(k) is 


( 10 ) 


E(k) = M B gH 0 + 2( S z ) J(0) jl - (11) 

MAGNETIZATION 

Application of a magnetic field is assumed sufficient to orient the net magnetization 
vector along the direction of the field. The field also introduces additional long-range 
order. The long-range order is thus due to the internal field and to the applied magne- 
tic field. The average z component of the spin is a measure of this long-range order. 
The thermal average ( S z ) is proportional to the magnetization. 

To calculate the magnetization ( S z ) , one uses the commutation relation of equa- 
tion (3) as well as 


S ■ S = S S + + S z + (s z ) 2 


( 12 ) 


which fixes the magnitude of the spin. The average value of the spin correlation is 
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(13) 


<S“S + > = S(S + 1) - <S Z > - (S z ) 


2 


For spin 1/2 equation (13) reduces to 


< S“S + > = - - < S Z > (14) 

2 

In this report only spin 1/2 is considered. With the use of equation (10) one can calcu- 
late the correlation function (S^S*) from the relation 


<S m ( t ’)S-(t)> = lim i 



[ A gm( E + ~ A gm^ E ~ i g)] eiE dE(t ~ t,} 

e E ( k )/3 _ i ' 


(15) 


Equations (10), (14), and (15) enable us to obtain an expression for magnetization. 

The calculation of the correlation function follows directly from equation (15) in the 
limit t = t\ From the 6 function representation 


6(x) = — —) (16) 

2ffi \x - i| x + i£/ 

it follows that 

<S“S + >=2<S z )cp (17) 


where 



(18) 


in which f} = l/(k B T) and E(k) is given by equation (11). Using equations (14) and (17), 
one finds that 


<S‘S + ) = 2<S z )cp =- - < S z ) 
2 


(19) 
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The expression for (S z ) then becomes 


1 

<S Z > 

1 + 2cp 


( 20 ) 


and it is easily shown that 


iV UjiM .il pi, 

2n4-(L 2 J 

k 

ANALYTICAL SOLUTION FOR MAGNETIZATION 

To calculate the magnetization as a function of temperature, the sums in equa- 
tion (21) must be evaluated over all values of k in the first Brillouin zone of the appro- 
priate lattice. Except at the very low and high temperature limits, numerical methods 
are usually used. However, such numerical solutions are difficult. One of the purposes 
of this report is to show that analytical solutions of equation (21) are possible if the 
methods of references 3 to 6 are used. 

Consider only crystals with cubic symmetry such as body-centered cubic (bcc). 

For such a lattice one can replace the sum which appears in equation (21) by an integral 



( 22 ) 


where for a bcc lattice 


x = 


k x a 




v 

y = - J — 

2 


z = 


V 

2 J 


(23) 
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where a is the length of the cell edge. Rewrite equation (20) in the form 


<s z > «-L 

2 A 


( 24 ) 


where 


_ 1 V coth /3E(k) 


A = — 
N 


4'’ 2 


(25) 


Using equation (22), one can write equation (25) as 


-/// 

n a «'0 «A) 


coth /3E(k) ( j x 


For the case of the bcc lattice, 


J(k) _ / cos V \ / cos V \ / cos k z a \ 


Using the identity 


(26) 


(27) 


coth nt = 


OO 

=i t ?y -j— 

irt t t Z— / R 2 . t 
R=1 K +l 


(28) 


and the techniques developed by Flax (ref. 5) and Flax and Raich (refs. 3 and 4), one can 
obtain 


A = — 

[2K(k)] 

2 2 r 

+ coth P - - + Q- csch 2 P coth P - 

/l\ 3 ‘ 

P 

- 7 T J 

P 8 

\P/ 


( 29 ) 


where 
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r 


p=® + Q=Q+® 
2 2 


( 30 ) 


^B gH 0 _ IT 
J(0)r t 


(31) 


Q = <^> 

r 

k R T 

T "j(0)' 

K(k) is a complete elliptic integral of the first kind, and 



(32) 

(33) 


(34) 


The solution of equation (29) gives A as a function of temperature, magnetic field, and 
magnetization. Hence, when equation (29) is substituted into equation (24), the magneti- 
zation is obtained. 

In the absence of an applied magnetic field the random phase approximation predicts 
that the magnetization shows a second-order transition which has a critical temperature 
with the value 


T — 0. 3589 J (0) _ q) 

k b 

It is then feasible to relate Hq and H' as follows: 


H’ = 


MggHo 

J(0) 


0.3588 M B gH Q 


k b T 


c 


(35) 


or 
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(36) 


h Q = 


k. T H' 
b c 


0.3588 




Substituting the appropriate values of the Boltzmann value, Lande g factor, and the 
Bohr magneton, one obtains 


1. 143X10 4 T H’ , 

H n = - 5. TlXlO" 5 T H’ 

u 2 c 


(37) 


Although the theory is for spin 1/2, and neither iron nor nickel has this spin, as exam- 
ples of the actual values of field corresponding to the nondimensional field H’, consider 
the following. Using the values for the critical temperature of iron (1043 K) and nickel 
(631 K), one can write 


H q = 596X10 4 H’ (iron) 

► 

H q = 365. 6xl0 4 H’ (nickel) 


(38) 


-4 -3 

The values used for H' are 5x10 and 5x10 , which correspond in the case of 

iron, to 2480 and 24 800 gauss, respectively. In the case of nickel these values corre- 
spond to 1802 and 18 020 gauss, respectively. 

Figure 1 shows a plot of ( S z ) as a function of temperature for several values of H' 
for a bcc lattice. When H' equals zero, there are no solutions for (S z ) above T . 
Thus, at this point the long-range order disappears, and above this point the present 
theory predicts complete disorder of the spins. When H' is not equal to zero, the fer- 
romagnetic transition occurs not at a single temperature but over a range of tempera- 
tures which forms a "Curie region, " so that the transition is broadened and results in 
the appearance of a "tail. " 

The reason for the tail is that the spins feel the effect of an applied magnetic field 
and the long-range order persists. The higher the field, the greater the broadening of 
the transition. 

Figure 2 shows a plot of (S z ) as a function of H' for several values of r for a bcc 
lattice. As can be seen from the figure, in the immediate neighborhood of the zero-field 
Curie point (r c = 0. 359) the magnetization has a nonlinear dependence on the field for all 
values of the field strength. When the temperature is raised, the nonlinearity remains 
in the weak-field region, but in the strong field there is a linear section of the magneti- 
zation curve. As one increases the temperature, the linear section of the curve extends 
farther toward weaker fields, and finally the magnetization curve goes entirely into a 
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straight line slope. The linear sections are simply Curie-law behavior. Experimental 
observations (ref. 12) show linear and nonlinear regions corresponding to those just 
described. 


MAGNETIC SUSCEPTIBILITY 


Let the magnetic susceptibility x be defined as 


x = (39) 

H’ a 

Using equations (24) to (34) in conjunction with equation (39), one obtains 


X = 


1 

2a A't 


(40) 


where 


A* =- 

"2K(k 1 r 

2 2 2 f 

+ coth P - — + — csch^P coth P - (- 

if] 

P 

- n 

P 8 L V 

p / _ 


(41) 


and 


P =-2! (1 + 2)0 
2 



(42) 


(43) 


For a given ferromagnetic material the spontaneous magnetization can occur only 
below a critical temperature T . Well above T c such materials are found to behave 
paramagnetically, and to have a well-defined susceptibility which follows the Curie- 
Weiss law, namely, 


X = 



(44) 
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in which C is the Curie constant and 6 is called the paramagnetic Curie temperature. 
Experimental measurements (ref. 12) show, however, that in the immediate neighbor- 
hood of the Curie point, the dependence of 1/x on r above T c is not linear. 

A plot of 1/x as a function of r from equation (40) is shown in figure 3. As re- 
quired the curve becomes linear at temperatures sufficiently far above the Curie point. 
The nonlinearity in the neighborhood of the Curie point agrees qualitatively with exper- 
iment (ref. 12). Previous theories such as molecular field theories have been unable 
to demonstrate this nonlinear behavior. 

CONCLUSIONS 

The thermodynamic properties of the Heisenberg ferromagnet in a magnetic field 
were investigated theoretically by using a random phase approximation. Expressions 
were derived from which the thermodynamic parameters can be calculated, and the 
following conclusions were drawn: 

1. The random phase approximation predicts that in the presence of a magnetic 
field there is a broadening of the ferromagnetic transition which results in the appear- 
ance of a tail in the magnetization-temperature curve. A higher magnetic field in- 
creases the broadening of the transition. 

2. There is a nonlinear dependence of the magnetization upon field at the Curie 
point. When the temperature is raised, the nonlinearity remains in the weak-field re- 
gion, but in the strong-field region the magnetization becomes linear with field, in 
agreement with the Curie law. 

3. Analytical expressions for magnetization and susceptibility are obtained which 
are valid for all temperatures. The analysis predicts that well above the Curie point 
the susceptibility follows the Curie-Weiss law. In the immediate neighborhood of the 
Curie point, the dependence of the reciprocal of susceptibility 1/x on temperature is 
nonlinear. The results are in qualitative agreement with experimental data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 7, 1971, 

129-02. 
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APPENDIX - SYMBOLS 


c 

Curie constant 

R 

sum index 

E k 

excitation energy 

S i 

spin operator at site i 

g 

Lande g factor 

<s z > 

magnetization 

H 

Hamiltonian 

T 

temperature 

H 0 

applied magnetic field 

T c 

Curie temperature 

H' 

M B gH 0 /J(0) 

a 

defined by eq. (31) 

i, i 

lattice sites 

P 

l/(kT) 

J(k) 

defined by eq. (9) 

T 

reduced temperature 

J i,J 

exchange interaction 

X 

<S Z )/H’ = Q/a 

K(k) 

elliptic integral of first kind 

<< )) 

Fourier transform of Green's 

k b 

Boltzmann's constant 


function 

P 

(a/2) + Q 

< ) 

conical average of operator 
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Magnetization, (s z ) 



Figure 1. - Magnetization as function of 
reduced temperature. 
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Figure 3. - Reciprocal of susceptibility as function of 
reduced temperature. 
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